Negative differential Rashba effect in two-dimensional hole systems 
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We demonstrate experimentally and theoretically that two-dimensional (2D) heavy hole systems 
in single heterostructures exhibit a decrease in spin-orbit interaction-induced spin splitting with an 
increase in perpendicular electric field. Using front and back gates, we measure the spin splitting 
as a function of applied electric field while keeping the density constant. Our results are in contrast 
to the more familiar case of 2D electrons where spin splitting increases with electric field. 
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In a solid that lacks inversion symmetry, the spin-orbit 
interaction leads to a lifting of the spin degeneracy of the 
energy bands, even in the absence of an applied magnetic 
field, B. In such a solid, the energy bands at finite wave 
vectors are split into two spin subbands with different 
energy surfaces, populations, and effective masses. The 
problem of inversion asymmetry-induced spin splitting 
in two-dimensional (2D) carrier systems in semiconduc- 
tor hetero junctions and quantum wells 0, |3, tJ |j| has 
become of renewed interest recently Q because of their 
possible use in realizing spintronic devices such as a spin 
field-effect transistor [flQl, an d for studying fundamental 
phenomena such as the spin Berry phase ||, ■ 

In 2D carrier systems confined to GaAs/AlGaAs het- 
erostructures, the bulk inversion asymmetry (BIA) of the 
zinc blende structure and the structure inversion asym- 
metry (SIA) of the confining potential contribute to the 
B = spin splitting [j, y| • While BIA is fixed, the so 
called Rashba spin splitting ^} due to SIA can be tuned 
by means of external gates that change the perpendic- 
ular electric field (E±) in the sample. For many years 
it has been assumed that the Rashba spin splitting in 
2D carrier systems is proportional to Ej_ that charac- 
terizes the inversion asymmetry of the confining poten- 
tial Q. 2D holes contained in a GaAs square quantum 
well provide an example On the contrary, in the 

present work we show both experimentally and theoret- 
ically that for heavy holes confined to a triangular well 
at the GaAs/AlGaAs interface, spin splitting decreases 
with an increase in Ej_. We demonstrate this negative 
differential Rashba effect by analyzing the Shubnikov-de 
Haas oscillations in this system at a constant density. 
We note that hole systems have recently gained great 
attention for spintronics applications |l2j | because ferro- 
magnetic (III,Mn)V compounds are intrinsically p type. 
A detailed understanding of the B — spin splitting in 
hole systems is thus of great importance. 

The sample used in our study was grown on a 
GaAs (001) substrate by molecular beam epitaxy and 
contains a modulation-doped 2D hole system confined 
to a GaAs/AlGaAs heterostructure [Fig. da)]. The 
Alo.3Gao.7As/GaAs interface is separated from a 16 nm 
thick Be-doped Alo.3Gao.7As layer (Be concentration of 
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FIG. 1: (a) Schematic sample cross section, (b) Schematic 
demonstrating how the gate voltages change the shape of the 
2D heterostructure potential (lines) and the charge density 
profile (shaded) while keeping the density constant. 



3.5 x 10 cm~ 3 ) by a 25 nm Alo.3Gao.7As spacer layer. 
We fabricated Hall bar samples via lithography and used 
In/Zn alloyed at 440°C for the ohmic contacts. Metal 
gates were deposited on the sample's front and back to 
control the 2D hole density (p) and E±. The low tem- 
perature mobility for the sample is 7.7 x 10 4 cm 2 /Vs at 
p = 2.3 x 10 11 cm -2 . We measured the longitudinal (R xx ) 
and transverse {R xy ) magneto-resistances at T m 30 mK 
via a standard low frequency lock-in technique. 

In single heterostructures, where SIA is the dominant 
source of spin splitting, the electric field E± experienced 
by the carriers is determined by the density-dependent 
self-consistent potential. This potential is determined, in 
turn, by the sample structure (spacer layer thickness and 
doping, etc.), and the applied gate biases. In our mea- 
surements, we used front and back gate biases to change 
the potential's profile and hence E± [Fig. Epa)], while 
keeping the density constant. We initially set the front- 
gate voltage (Vfg) to 0.55 V and back-gate voltage (Vbg) 
to —100 V with respect to the 2D hole system, leading to 
p = 1.84 x 10 11 cm" 2 , and measured R xy and R xx as a 
function of B between -3 T to 5 T. Then at B = 1 T, we 
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FIG. 2: Observed Shubriikov-de Haas oscillations for a 2D 
hole system confined to a (001) GaAs/AlGaAs single het- 
erostructure at two different Ex- Inset: The Fourier spectra 
of these oscillations at the corresponding electric fields. The 
dotted curves in the main figure and the inset are shifted ver- 
tically for clarity. 



FIG. 3: Spin splitting, A/, versus the change AE± in the 
perpendicular electric field for the measured (o = /tot — 2/_ , 
A = /+ — /_) and calculated (solid line) magneto-oscillations. 
The effect of varying E± on the shape of the 2D heterostruc- 
ture potential (lines) and the charge density profile (shaded) 
is also shown. Inset: The increase in the energy gap Af' 
between the first HH and LH subbands is shown. 



increased Vbg and noted the change in R xy ; this change 
in R xy gives the corresponding change in density (Ap). 
Vfg w &s then decreased to recover the original R xy and 
hence the original p. This procedure leads to a change 
AEx = eAp/e in E±_ (e is the electron charge and e is 
the dielectric constant) while keeping the density con- 
stant to within 1%. AE± is measured with respect to 
E x at V FG = 0.55 V and Vbg = -100 V. 

Figure shows the low-field Shubnikov-de Haas (SdH) 
oscillations for two Ex differing by 8 kV/cm. The Fourier 
transform (FT) spectra of these oscillations, shown in 
Fig- [21 (inset), exhibit four dominant peaks at frequen- 
cies /_, / avg , /+, and /tot, with the relation / tot = 
/+ + /- = 2/ aV g. The /t t frequency, when multiplied by 
e/h, matches well the total 2D hole density deduced from 
the Hall resistance (h is the Planck's constant). The two 
peaks at /_ and /+ correspond to the holes in individual 
spin subbands although their positions times e /h do not 
exactly give the spin subband densities jElUlil. Nev- 
ertheless, as discussed below, this discrepancy between 
(e/h) f± and the B = spin subband densities is minor 
and A/ = /+ — /_ = /tot — 2/_ provides a good measure 
of the spin splitting. The vertical lines in the inset of 
Fig. [3 at the /_ and /+ peaks clearly indicate that A/ 
decreases when AEx is increased from to 8 kV/cm. 
The vertical line at the /tot peak shows that the total 
hole density is held constant. 

We compare the experimental data with accurate nu- 
merical calculations of the magneto-oscillations at B > 0. 
First we perform fully self-consistent calculations of the 
subband structure at B = in order to obtain the 



Hartree potential Vh as a function of Ex 15] . We 
assumed that the concentration of unintentional back- 
ground impurities in the GaAs space charge layer was 
1 x 10 14 cm~ 3 17]. This assumption is based on our 
sample parameters; we note, however, that the deduced 
spin splitting is insensitive to the exact value of the back- 
ground doping. Using Vh we obtain the Landau fan chart 
for B > from an 8 x 8 k • p Hamiltonian that fully takes 
into account the spin-orbit coupling due to both SIA and 
BIA Ulig. From this fan chart we then determine the 
magneto-oscillations by evaluating the density of states 
at the Fermi energy as a function of B. 

In Fig.|3|the experimental and calculated A/ from the 
corresponding FT spectra are plotted versus AEx- It is 
clear that increasing Ex lowers the spin splitting [l8j . 
We also verified that the calculated B = spin splitting, 
defined as the difference between the spin subband den- 
sities, shows the same negative differential Rashba effect 
as A/ obtained from the calculated magneto-oscillations. 
The difference between the B = spin splitting times h/e 
and A/ is only < 0.14 T in the range of Ex shown in Fig. 
13 This confirms that the B — spin splitting shows the 
same negative differential trend. 

We can understand these surprising results in the fol- 
lowing way. The hole states in the uppermost valence 
band Tg have the angular momentum j = 3/2. In 
2D systems, the four hole states split into heavy-hole 
(HH) states with z component of angular momentum 
to = ±3/2 and light-hole (LH) states with to = ±1/2. 
Here the quantization axis is perpendicular to the 2D 
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plane. On the other hand, the Rashba spin-orbit coupling 
acts like a fc-dependent effective magnetic field which is 
oriented in the plane so that it favors to orient the quan- 
tization axis of the angular momentum in-plane. How- 
ever, this is not possible within the subspace of HH states 
(m = ±3/2) so that — in contrast to j — 1/2 electron 
systems — the Rashba spin splitting of HH states is a 
higher-order effect. Neglecting anisotropic corrections, it 
is characterized by the Hamiltonian p| 
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with a± — 1/2((T X ± i<Jy) and k± = k x ± ik y , where a x 
and o v are the Pauli spin matrices in the x and y direc- 
tions respectively. Using third-order Lowdin perturba- 
tion theory we obtain for the Rashba coefficient [5^ of 
the lowest HH subband a = 1 

"Tl -A | i V -i, » -A,o / -A,o -A | 

(2) 

where 72 and 73 are the Luttinger parameters [16j and 
A^, = £ v a — £ v a , where and £ l a are the energies of 
the ath HH and LH subband, respectively. The symbol 
a denotes a numerical prefactor which depends on the 
geometry of the quasi-2D system. We can estimate the 
value of a assuming an infinitely deep rectangular quan- 
tum well which yields a — 64/(97r 2 ). We see from Eq. 
(0) that the Rashba spin splitting of the HH states de- 
pends not only on the electric field E± but also on the 
separation between the HH and LH subbands. As can 



be seen in Fig. [3] (inset), the gap between the HH-LH 
subbands increases with an increase in E± , giving rise to 
a decreasing Rashba coefficient This result reflects 
the fact that a large HH-LH splitting yields a "rigid" an- 
gular momentum perpendicular to the 2D plane so that 
the Rashba spin splitting will be suppressed. 

We can estimate the effect of changing E± using the 
well-known triangular potential approximation jl9( . Here 
we have, for the subband energies ££ measured from the 
valence band edge, ££ cx E± 3 which implies oc E± 3 . 
Therefore, we can expect from Eqs. and J5J that the 
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Rashba spin splitting decreases proportional to E ± 
when E± is increased, in agreement with our more accu- 
rate numerical calculations. 

In a previous study [2(| this surprising behavior of 
PiE± was shown as a function of density. By lower- 
ing the density, spin splitting and E± both decrease but 
the term PiE± increases. However in the present work, 
by keeping the density constant and only varying E± , we 
are able to directly demonstrate the negative differential 
Rashba effect in heavy hole 2D systems confined to single 
GaAs/AlGaAs heterostructures. 

In conclusion, our study highlights the subtle and un- 
expected dependence of the Rashba spin splitting on E± 
in 2D hole systems. The results are important for the 
spintronic devices 0, 0, ^3 whose operation relies on the 
tuning of the spin splitting via applied electric field. 
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